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ABSTRACT 

We build a bulk velocity-dependent photoionization model of the warm ab- 
sorber of the Seyfert 1 galaxy NGC 3783. By adopting functional forms for 
velocity of the flow and its particle density with radius, appropriate for radiation 
driven winds, we compute the ionization, temperature, line Doppler shift, and 
line optical depths as a function of distance. By doing this we obtain detailed pro- 
files for the entire absorption line spectrum. The model reproduces the observed 
relationship between the gas ionization and the velocity shift of the line centroids 
as well as the asymmetry of the absorption lines in the X-ray spectrum of NGC 
3783. It is found that the distribution of asymmetry requires the presence of two 
outfiows: a higher ionization component responsible for the blue wings of the 
high ionization lines and the red wings of the low ionization oxygen lines, and a 
lower ionization fiow that makes up for the blue wings of the oxygen lines. Our 
model predicts a relationship between the X-rays and the UV absorbers, where 
the component creating the long wavelength lines in the X-ray spectrum is also 
responsible for creating at least part of the UV absorption troughs. 

Subject headings: AGNs - Photoionization - Plasma - Asymmetry - Absorption 
lines 



Introduction 



X-ray (0.5-10 keV) Chandra spectra of type-I active galactic nuclei (AGN) often show 
presence of absorption lines coming from H- and He-like ions of O, Ne, Mg, Si, S as well 
as from Fe xvil-Fe xxill L-shell transitions. These lines are near the region of the bound- 
free absor ption edges of O vn and O Vlll at ~ 0.8 keV, which are the hallmark of warm 
absorbers (IGeorge et al.lll998l : lKomossalll999l ). In Seyfert 1 galaxies these absorption lines 
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can be narro w (Seyfert 1 NAL s) with full width at half maximum (FWHM) spanning ~ 100- 
500 km s"^ kaspi et al.l boooh. or broad (Seyfert 1 BALs) with FWHM ~ 500-2000 km s"^ 



( iKaspi et al. 



The troughs of the lines are blue-shifted relative to rest fra me of the host 



galax y by velocities that go from several hundred to a few thousand km s ^ (lYaqoob et al. 



20031 ) , i mplying the presence of outflows that cover a wide range of velocities and ionizati on 
stages jKaastra et all I2OO0I : iKaspi et all I2OO2I : iKrolik fc Kris] I2OO1I : iRozahska eraP booeh . 



In Seyfert galaxies and Quasars absorption spectra show a range of ioniz ation stages 
that arise from ionization parameters that span ^ ~ 10 — 1000 erg cm s~^ (IKaspi et al. 
20021). Simultaneous UV and X-ray obser vations may enhance this range even further 



([Shields fc HamannI 119971 : iGabel et al.ll2005l ). The location of these absorbing material is 
uncertain. Models of X-ray absorbers in AGN place them in a wi de range of distance s from 
the c entral source, from winds originate d at the accret i on di sk (IMurray et al.l Il995l : lElvis 
2OOOI ) ■ out to the dusty (~ 1 pc) torus (IKrolik fc Krissll200ll ) and beyond the narrow-line 
region (e.g. lOgle et al.ll2000l ). 



In the Seyfert galaxy NGC 3783 the 900 ks Chandra spectrum of iKaspi et al.l (120021 ) 
allows for precise measurements of radial velocities and widths of the lines. It is seen that 
the velocity shift of lines from Fe xxill-Mg xil cover a range of ~ 60 — 600 km s~^ while the 
lowly ionized lines Si xill-0 vil cover velocities ~ 500 — 1000 km s~^ (see Figure 6 in Ramirez 
et. al. 2005). The average velocity of the warm absorber outflow of NGC 37 83 is around 
~ 500 km s~^. The spectrum also reveals that the line profiles are asymmetric (IKaspi et al. 



2OO2I ) . in such a way that approximately 90% of the lines have extended blue wings. Such 
asymmetries were quantified by Ramirez et al. (2005). In terms of ionization, mostly high 
ionization species are seen in the short-wavelength portion of the spectrum ~ 4 — 12 A. 
Here, resonant hues from Fe xxill, Fe xxil, Fe xxi, S xvi, S XV, Si xiv, and Mg xil cover 
ionization parameters ^ from ~ 630 to ~ 150 ergs cm s~^. The longer wavelength of the 
spectrum, ~ 10 — 20 A, is dominated by lower ionization lines from Si xill, Fe xvill, Ne x, 
Mg XI, Fe XVII, Ne ix, O vill, and O vil that span ionization parameters ^ from ~ 125 ergs 



cm s 



down to 



ergs cm s 



Thus, photoionization modeling shows that the observed spectrum is hardly explained by 
a single ionization parameter. Rath er two or three components are needed, and this has been 



the s ubject of well-detailed studies ( iKrongold et al.ll2003l : iNetzer et al.ll2003t iKrongold et al. 



2OO5I ). In those works two/three components in pressure equilibrium are enough to account 
for all the charge states seen in the spectra. When start runing those simulations, the ion- 
ization parameter is not constant throughout the cloud but it varies in a self consistent way. 
On the other hand ISteenbrugge et al.l ( 120031 ) used a nearly continuous distribution of ion- 
ization parameters that span three orders of magnitude. In reality both approaches do not 
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need to be necessarily contradictory. We just have to look at t he po ss ibility of intr o ducin g 
species between the c omponents proposed by iKrongold et al.l (120031 ): iNetzer et al.l ( 120031 ): 
Krongold et al.l (|2005[ ): and that immediately lead us to the discussion of whether the ma- 
terial absorbing is clumpy or not. And that is in part the goal of the present work. We 
present both possibilities, continuous vs cumply flows, with their physical implications along 
the article. 

Although there is no evide nt correlation between the ionization of t he absorption species 
and the velocity shifts of lines (IKaspi et al.ll2002l ). iRamirez et al.l (120051 ) open the possibility 
for such relationship to exist, as the lines are consistent with being originated from a n 
outward flowing wind. In that work a wind velocity law was adopted (ICastor et al.lll975l ). 
with an ionization fraction that goes as a power-law of the ionization parameter q ^ with 
1] > for all C, and all l ines. Further, q ^ and an optical depth of the form proposed 



by the SEI method from lLamers et al.l (119871 ) were used to fit the line in the spectrum and 
study the possible relationship between the profiles and the expansion velocity of the flow. 

Here, in this theoretical approach no analytical dependence of q with ^ is assumed, but 
computed self-consistently using the photoionization code XSTAR with modifications of the 
optical depth of the lines for working with expanding outflow. 

It is the aim of the present paper to further study the wind scenario for the warm 
absorber of NGC 3783 by modeling a photoionized wind flow and trying to reproduce the 
main features in the spectrum in terms of both equivalent widths and line profiles. The model 
consists of a spherically syr nmetric gas fl o w wit h increasing velocity with radius according 
to the wind velocity law of ICastor et al.l ( 119751 ). The microphysics of the plasma is then 
solved in detail along the wind, and radiative transfer is treated for the flowing plasma 
including Doppler shift effects on the emerging spectrum. This kind of approach has been 
used p reviously in studies of t he ionization and thermal properties of O stars by iDrew 
(1989), cataclysmic variables by Drew &: VerbuntI (1985) and in an evaluation of absorption 
line profiles from winds in AGN by iDrew &: Giddingsl ( 119821 ). 



e.g 

(1- 



Self-consisterit hyd rodynamic modeling of AGNs have been performed in the past (see 



Murray et al.lll995l ) which qualitatively predict the same /3- velocity law, where v{r) 



tq/tY , as adopted here. However, these models did not include detailed treatment of 
thermal and spectral processes that lead to synthetic spectra. 

The present paper is organized as follows: In §5] we discuss the statistical significance 
of the X-ray lines under study. Later in §21 we present the theoretical method we have used 
in this work. In f|5] we present the main results of our work. We present the final physical 
solution in ^ And we summarize in ^ 
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2. Significance of the shifting in the X-ray Lines 



Since our analysis searches for possible correlation between ionization state of the ions 
forming the absorption lines observed in the X-ray spectrum of this object, and the Doppler 
velocity shift of those lines, we first of all study the statistical significance of this possible 
shifting. A convenient way to do that is through the separation of the ions into groups and 
look for th e velocity shi f ting. We take all the ions we find in Table [31 The line's centroid are 
take n fromlKaspi et al.l (120021). T he ionization fraction curve s for the classification are taken 
from lKallman fc Bautistal ( 1200 ll ) ; iBautista fc Kallmaru (120011 ) . Afterwards we grouping them 
into three groups 0: Group (1), representative of all ions of low ionization with ionization 
parameter < ^(ergs cm s^^) < 10; Group (2) 10 < ^(ergs cm s~^) < 250 representative of 
intermediate ionization state; and Group (3) 250 < .^(ergs cm s~^) < 650. A graphical way 
to represent this grouping scheme is with a boxplot shown in Fig [H Each box is made of 
five-number summaries: the smallest observation (sample minimum, which are the extreme 
thin bars at the left of each box), a lower quartile (Ql, which is the left thick border of each 
box), median (Q2, black thick vertical line), upper quartile (Q3, the right thick border of 
each box) and largest observation (sample maximum, which are the extreme thin bars at the 
right of each box). The first conclusion is that all the three Q2s are different. It is clear that 
Groups 1 and 3 are different at > Icr. There is overlap between Groups 1 and 2 (possible 
due to limited resolution of the instrument), but Groups 2 and 3 are different at > Icr. We 
go ahead in our analysis bearing in mind these overlaps in groups. 



3. Method 



Let us consider a radiation source with L 
massive (~ IO^Mq) black hole (BH). Material ^ 



, 10^^ — 10^^ ergs s ^, arising from a super- 
0.1-1 pc from the BH only needs to absorb 



(Arav et al. 


1994; 


Ramirez 


2008; 


Saez et al. 



20091 : IChartas et al.ll2009[ ). By conservation of mass the number density of hydrogen can be 



written in spherical symmetry as 



M 



(1) 



where M is the mass-loss rate, v{r) is the outflow speed at radius r, fi is the mean atomic 
weight per hydrogen atom and mn is the hydrogen mass. 



^We group and generate the figure using the statistical package R - 



http : / / www . r-pro j ect . org 
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We adop ted a velocity law v (r) compatible with the predictions of the radiatively driven 
wind theory (jCastor et al.lll975[ ). The velocity law has two fundamental roles. The first is 
to shift the frequency of the absorbing lines according to the Doppler effect. The second 
is to dilute the gas density, affecting radiative transfer across the gas and consequently the 
ionization and thermal state of the gas. The velocity law varies with distance as 



w{x) = Wo + {l-Wo)[l 



X J 



(2) 



here w is the velocity normalized to the terminal velocity of the wind foo, is the velocity 
in the base of the wind and x is the distance normalized to the radius of the central core 
r/rg. The parameter /3 is the quantity governing the slope of the velocity with the distance 
and its ad hoc value depends on the ty pe of radiative force acting on the wind. Analysis 
of hot stars suggests th at 0.5 5, /3 5, 1 (ILamers et al.l 119871 ). The evaluation of radiation- 
drive n wind in AGNs of iDrew fc Giddings (1l982l ). and more recent dynamical calculations 



[e.g, 



Murray et al.lll995l : iProga et al. 



20001 ). suggest that in AGNs 0.5 < /3 < 2. We have 
computed models using velocity laws /3 = 0.5, 1, 1.5 and 2, representing from fast (/3 = 0.5) 
to slow winds (/3 = 2). We found that fast winds were not able to simultaneously cover the 
range in velocity and ionization state observed in the spectrum of NGC 3783, needed for the 
goal of this work, and that the best-fit velocity law was /3 = 2. This is why for the rest of 
this work we use /3 = 2. Then, we rewrite the number density in terms of the velocity law 
as, 

nnix) = tiqx^'^w^'^, (3) 



where no = [j^^Tq \ 
Doppler relation. 



where Aq is rest wavelength. 



The absorbing line frequencies are shifted according to the 



w 



c 



1 - 



A 



(4) 



Our models are based on clouds illuminated by a point-like X-ray source. The input 
parameters are the source spectrum, the gas composition, the gas density nnix), and the 
outflow velocity w{x), where x is the position of each slab inside the cloud normalized to the 
radius of the most exposed face to the source, tq. The source spectrum is described by the 
spectral luminosity = Lf^, where L is the integrated luminosity from 1 to 1000 Ryd, and 
j^ioooRyd j.^^^ _ ^]^[^ spectral function is taken to be a power law ~ e^", with a = 1. 
The gas consists of the foll owing elements, H, He , C, N, O, Ne, Mg, Si, S, Ar, Ca and Fe. 
We use solar abundances of iGrevesse et al.l (119961 ) . in all our models. 



T hermal and statistical equilibrium in our models ar e computed with the code XS- 
TAR (jKallman &: Bautista 2001; Bautista fc Kallman 2001). The code includes all relevant 
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atomic processes and computes the equilibrium temperature and optical depths of the most 
prominent X-ray and UV lines identified in AGN spectra. 

We consider two types of models, the single absorber model (SA) and the multicompo- 
nent model (MC). In both cases the absorption fine profiles depend simultaneously upon the 
ionization and the kinematics of the absorbing gas. 



3.1. The Single absorber Model 



This model consists of a single extended cloud directly in the line of sight between the 
observer and the central source. Such a cloud fiows away from the central source and towards 
the observed. For t his model we make use of the Sobolev approximation for computing the 
line optical depths (ICastor et al.l 119751 ) . 



7re 



T,Ar] 



f dv\ ^ 



mc \dr J 



(5) 



where / is the absorption oscillator strength, Aq (in cm) is the laboratory wavelength of the 
transition, rii (in cm~^) is the number density of the absorbing ion, and dv/dr is the velocity 
gradient in the wind. This gives us the relation between the outfiow state and the radiation 
field. This is different from the calculation of the optical depth in the static case, which is 
directly proportional to the column density. 

Inside the cloud we us e a one-step forward differencing formula for the radiation transfer 
kallman fc McCravlE982l ) 



where ju{r) is the emission coefficient at radius r. 



1 _ e-rAr) 



(6) 



Under ionization equilibrium condition the state of the gas depends just up on the shape 
of the ionizing spectrum and the ionization parameter ^, that we define as in iTarter et al. 
fll969[ ) 

47rF(r) 



nnir) 



(7) 



where F{r) is the total ionizing fiux. 



F{r) 



lOOORyd 



Ly{r)dv, 
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with -Ljy(r) given by equation ([6]). The requirement that C, spans various orders of magnitude, 
as observed, yields a cloud geometrically thick throughout most of the spectrum, i.e. a 
photoionization bounded cloud. For instance, a luminosity L ~ 10^^ ergs s~^and A,^ = 
.^2—^1 = -(^ — t) = 10'^ yield n ~ 10^ cm~^ and a column density of the absorbing 
material is Nh ~ nAR ~ 9.7 x 10^^ cm~^. This is a large value, such that even lines with 
moderately small oscillator strengths become saturated in the emergent spectrum, unless the 
metal abundances are reduced by several orders of magnitude with respect to solar. 

In Figure [2] we compare the computed optical depth for the Ne x A 12. 134 line in a 
stationary nebula with L = 10^^ ergs s~^, and constant density n = 10^ cm~'^ with a flow 
with n{x) = 3.3 x 10^x~^'^w~^, a wind velocity function with /3= 2, and Voo = 1000 km s~^. 
In both cases we adopt a turbulence velocity of 200 km s~^ and solar abundances. In the 
stationary case, r is proportional the column density of the absorbing material and can reach 
very large values. In the outflow model, r peaks at log^ ~ 1.3 {v ~ 890 km s~^) reaching 
~ 5000, that is nearly two orders of magnitude smaller than the stationary cloud for the same 
^. Such high optical depths are common to other lines in the X-ray band, like O vill in the 
14-20 A and the line O vill A18.969. In s uch cases lines appear saturated in the absorption 
spectrum in contrast with observation. In lRamirez et al.l (120051 ) we fit the integrated optical 
depths of resonant lines in the spectrum of NGC 3783 and found Ttot ^ 1 for most of the 
lines. We illustrate in Figure [3] the consequences of taking this model for the reproduction 
of the profile of the line Ne x A12.134. In order to have a non-saturated line (as observed), 
we had to reduce the abundance to ~ 1 % solar, which has no physical motivation. We do 
not go further with this model, and present the clumpy (multicomponent) scenario in the 
next section. 



3.2. The Multicomponent Model 



Now, we examine the scenario in which the absorption profile seen in the X-ray and UV 
spectra of AGNs are made up of multiple components. The main differerice between this 
model calculations and prev ious ones by other authors (e.g., Kaspi et"alll2002 : Krongold et al. 
20031 : iRozahska et al.ll2006l ). is that in our model all the components are linked by a velocity 
law and a gas density distribution. 

Each absorber is specified by an ionization parameter ^, a column density A'^^^, an 
absorption covering factor, a gas density nnir) and an outflow velocity v{r) in order to shift 
the absorption lines according to the Doppler effect. Once the ionic fraction is calculated 
from the ionization equations, and the ionic levels computed, the opacity in each frequency 
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bin is 

nugi 



f^i^P = fini 



0(Az/), (9) 

rrieC I riiguj 

where k^, is the opacity at the frequency z/, p = fimnnH is the mass density, fi is the 
absorption oscillator strength, ni, Uu (in cm~^), gi, gu are the number density and the 
statistical weight of the lower and upper levels of the transition respectively. We allow for 
the lines to have a finite width characterized by the line profile 0(Ai/), with a width which is 
the greater between the thermal and the turbulent motions. In all our models the turbulence 
velocity is assumed to be 200 km s~^. The optical depth of a line in each component is 

Ki,{r)p{r)dr, (10) 



where r2 and ri are the limits of the cloud. 

It is important to highlight a special difference between the way we use the ionization 
parameter in the MC model and the way it is used in the SA model. Because of the clouds 
intervening in this model are optically thin, the ionization parameter at the most exposed 
face of the cloud remains essentially constant through the cloud, i.e. 

e(r) = ^, (11) 

nH\r) 

where Fir) is the total ionizing flux at the radius r. 

When two or more outflows are put together they are assumed to be distributed such 
that the observed spectrum is the result of the addition of all components. So the radiative 
flux in each bin of frequency for the composite spectrum is 



where m is the number of absorbing clouds in the line-of-sight, fu{xi) is the flux resulting 
from the pass of the continuum radiation through the absorbing cloud "i", and Xi is the 
normalized spatial radius of the cloud. 



In iRamirez et al.l (120051 ) is showed that in order to fit the line profiles in NGC 3783 a 
geometry different from the spherical for the gas distribution is required. Such deviation 
from spherical geometry has two effects. In the normalized notation the number density is 



Unix) = nox-^+^w-\ (13) 

where < k < 2. A positive value of n implies that the gas flow dilutes more slowly than in 
a free spherical expansion, i.e. that there are sources of gas embedded in the flow, or that 



- 9 - 



the flow is confined. A negative value corresponds to sinks of gas in the flow, or expansion 
of an initially confined flow in a flaring geometry. Secondly, we allow the radiation flux to 
have a form 

F = FoX x-^-P, (14) 

where p is an index to mimic a deviation of flux from the pure geometrical dilution case. 

This is expected if the medium between clouds has a significant optical depth (if p is 
positive) or if there are sources of radiation embedded in the flow (if p is negative). 



4. Assumption about the number density 



Before going into the presentation of the results we would like to highlight important 
differences between the underlying physical/geometrical motivation of the present work and 
previous ones. After reviewing single (geometrically thick model) and clumpy (geometrically 
thin model), we favored the latest one and require for the production of t he line profile found 
in the spectrum of Seyfert galaxies, AR/R « 1 (also based on results by lGabel et al.l (120051 ) 
from UV data). 



Based on UV CIII* density constraints (IGabel et al.ll2005l ). the electron density for the 
absorber could be Ue ~ 10"^ cm~^. With this (using eq. [TT] and uh ~ ne/1.2), the distance 



betwe en the absorber and the central source is i? ~ 25 pc (^ 10 cm). As in iGabel et al. 



(120051 ) . the low ionization species (XLI in that paper) seems to share the same kinematics 
with the UV absorber. Using typical luminosity L for this object of ~ 10^^ ergs s~^, uh ^ 10^ 



cm ^, ionization parameter of ^ ~ 1, i? ~ 10^*^ cm; similar to that computed by IGabel et al. 
(120051 ) . but if we take uh ~ 10^^ cm~^, as is r equire d for AR « R; then R 10^^'^ cm. This 



is not in contradiction with iKrongold et al.l (120051 ). where they set limits on the density and 
location of the absorber; with > 10^ cm~ ^ , and D < 5.7 pc. However is clearly different 
from the computation made by lNetzer et al.l (120031 ). of ng < 5 x 10^ cm~^ D > 3.2 pc for the 
logiUox) = —2.4 component, Ue < 10^ cm~^ D > 0.63 pc for the logiJJox) = —1.2 component, 
and rig < 2.5x10^ cm~^D > 0.18 pc for the /o(y'(t/oa;) = —0.6 component. It is clear that some 
of the differences in the est imations can be due t o differenc es on the nature of the physics 
behind the estimations. In iKrongold et al.l (120051 ) (and also iReeves et al.l (120041 ) for the Fe 
K shell), the fundamental assumption behind the computation of is that the absorber 



respon ds instantaneously to changes of flux in the ionizing source, while in iNetzer et al. 



(120031 ) the estimations are based on temperature derived models, average recombination rates 
and no response to continuum variations on timescale of 10 days. And, for our purposes, that 
basically translate in differences between the two proposed physical mechanisms; thermally 
accelerated winds, which from grounds has to consider sublimination radius for the material 
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not being evaporated, and radiatively accelerated wind with origin possibly in the accretion 
disk with subparsec scales; both competitor theories in the explanation of the origen of the 
warm-absorber outflows. 



Results and Discussions 



In Figure m we show the variation of the different variables governing the kinematics and 
the ionization structure of one of our models (model A). This is a wind with a velocity slope 
P = 2, K = 1.5 and p = 1.5. The number of kinematic component m is 11. The assumed 
abundance is solar as is given in Table [H and each component has a column density of = 
5 X 10^° cm~^. In this model the variation in the ionization parameter is log ^ = 3.5 — (—0.65) 
[erg cm/s], where the "launching" ionization parameter is defined as log^o = 3.5 [erg cm/s], 
and the variation in density lognn = 11-4 — 10.26 [cm~'^]. We can see the variation with 
distance of the input parameters w, flux, uh and ^, for model A in Figure IH In Figure O 
we plot the velocity shifts taken fror n the maxini u m liri e optical depths Tmax of our model 
versus the hne centroids measured by iKaspi et al.l (120021 ) for NGC 3783 (see Table 3 in that 
paper). 

We fit a linear model by robust regression using the M estimator (Huber 1964). The 
weights are included as the inverse of the variances, with the variance equal to Av = {vio + 
Vhi)!^,, and fjo, and Vhi are the differences betwe en centroid a nd lo wer and upper velocity 
limits of the measured lines given in Table 3 of iKaspi et al.l (120021 ). The best-fit slope is 
0.97±0.31 (solid line in the figure), and the residual standard error (rms) is 15.4 for 25 
degrees of freedom (dof). This agreement is encouraging considering that no model have 
been suggested before to explain the possible correlation between the ionization parameter 
and the velocity shift seen in this Seyfert galaxy. 

Although the velocity shifts measured with respect to line minima are well explained by 
this model the distribution (in g eneral) of line profile shapes is not in agreement with that 
reported in iRamirez et al.l (120051 ) . for instance. From that study most lines have extended 
blue wings. On the other hand, the model does predict the correct asymmetries for a few 
lines from highly ionized species, but the profiles of the lower ionizations lines exhibit more 
extended red wings, unlike observations. The discrepancy between the modeled profile of 
the O VIII A18.9 line and the observed in NGC 3783, is important. This is why we created 
model B. 



In Figure [6] we present the variation of the variables governing the kinematics and the 
ionization structure of another model (model B). Here we change slightly the parameters 
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log^o, and foo (see Table [2] for details). The wind has the parameters (3 = 2, k = 1.5 and 
p = 1.5. The number of kinematic component m is 11. The abundances and the column 
density are as in model A. In this model the variations in ionization parameter and density 
are log^ = 3.00 — (—1.13) [erg cm/s] and logrin = 11.40 — 10.26 [cm~'^]. We compare 
models A and B in terms of the variation of optical depth vs. velocity of the flows of three 
important lines, i.e, Si xiv A6.182, Si xill A6.648 and Mg xil A7.106. We could estimate 
the contribution of each cloud to the formation of the composed profiles of these lines. We 
have that model B is superior to model A in describing the low-velocity portion of the flow 
reflected by the high ionization lines seen in NGC 3783 (extended blue wings), but yields a 
worse description of t he spectrum at hig her terminal velocities responsible for the center of 
the lines measured by lKaspi et al.l (120021 ). Figure [7] shows the fit of a linear model by robust 
regression using the M estimator for model B. As in model A, the weights are included as 
the inverse of the variances, with the variance equal to Av = {yio + Vhi)/^., and vio, and Vhi 
are the differences betwe en centroid a nd lo wer and upper velocity limits of the measured 
lines given in Table 3 of iKaspi et al.l ( 12002| ). The best-fit slope is 0.72±0.20 (solid line in 
the figure), and the residual standard error (rms) is 17.6 for 25 degrees of freedom (dof). In 
general, the scatter get worse for lines with terminal velocities greater than 1100 km s^^. 



5.1. Modeling global properties - ^ vs velocity 



One of the prime goals of this study is to examine the relationship between the kine- 
matics and the ionization structure of the flow. This is in com plement with earlier works 
(IKrongold et al.l l2003t iNetzer et al.l 120031 : iKrongold et al.l |2005|) which follow the algorithm 
of constructing a grid of static photoionization models varying in ionization parameter and 
column density with the selection of this quantities which best reproduce the EW observed 
in the spectra, in addition paying attention to the relationship between these physical pa- 
rameters and the position in wavelength space of the absorption troughs. 



Krongold et al.l (120031 ) fitted photoionization models to the 900 ks spectrum of NGC 
3783 using two phases; one at high ionization and high temperature logT ~ 5.98 (HIP) and 
one at low ionization and temperature logT ~ 4.41 (LIP), finding good agreement between 
the modeled equivalent wid th (EW) and the measured for a set of absorption lines (see Figure 
9 of IKrongold et al.l (120031 )). They set each of these phases at a single outfiow velocity of 



750 km s ^, assuming spati al coexistence of the two absorbers. The approach is similar to 



that from Netzer et al 



(120031 ) ■ but in the latest one they use three components in pressure 



equilibrium, with two kinematic components each. 
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One property of our model is the capability of open the possibility of establishing □ 
a relationship between the ionization and the kinematics of the gas based on a radiatively 
accelerated wind. Figures |H] and show this relation for models A and B respectively. In these 
plots show the predicted relationship between observed velocities and ionization parameter. 
While absorption from highly ionized ions originates from low-to-intermediate (200-600 km 
s~^) velocities, lines from lower ionization stages are formed at intermediate-to-high velocities 
(600-1000 km s'^). 



5.2. Far-UV and UV absorbers 



The relationship between the X-ray abs orption spectrum and the UV absorption spec- 
trurn is at present a subject of controversy. iKraemer et al.l (1200 ll ) and iGabel et al.l ( l2003l . 
20051 ) have analyzed de UV spectrum of NGC 3783. HST/STIS and FUSE spectra show ab- 
sorption troughs from the low order Lyman series (i.e. Lya, Ly/3, Ly7), C IV AA1548.2, 1550.8, 
N V AA1238.8, 1242.8, O vi AA1032, 1038. All these lines are seen in three kinematic com- 
ponents at —1365, —548, and —7 24 km s~^ (compon e nts 1 ,2 and 3 respectively). A weak 
fourth component is reported by iGabel et al.l ( l2003l . |2005| ) at —1027 km s~^. Figure [TOl 
shows the spectrum predicted by model B in the range 700-1700 A. It is clear from this 
figure that lines due to low order Lyman series. He ll, C iv, N v, O VI, and Ne Vlll would 
be detectable in the UV band of the spectrum. Figure [TT] shows the spectra of several UV 
lines as a function of radial velocity respect to the systemic. At the top of the Figure we 
present the Lya blended with the He ll A1215 line (upper left) which yields a feature centered 
at ~ —1300 km s^^. Also an absorption feature is formed with the Ly/3 and He ll A1025 
lines (upper center) with center at ~ —1200 km s~^, and the line Ly7 (upper right) with 
a velocity ~ —1100 km s^^. At the bottom of the Figure we plot the spectrum of three 
important doublets (C IV AA1548.2, 1550.8, N v AA1238.8, 1242.8, and O vi AA1032, 1038). 
The solid line depicts the velocity spectra constructed taking the shorter wavelength of the 
doublet. The spectra in dashed lines were made taking the longer wavelength. One can 
see the similarity between the velocities predicted by our model and the high-velocity com- 
ponents (1, 4 and likely 3) seen in the UV spectrum of NGC 3783 and othe rs Seyfert 1 



galaxi es (for exar nple NGC 5548). T his is similar to the con clusion reached by iGabel et al. 



(120051 ) , based on iKaspi et al.l ( 12002| ) and iGabel et al.l (120031 ) , where all X-ray lines having 
sufficiently high-resolution and S/N were found to span the radial velocities of the three UV 
kinematic components (see also next section). These are predictions which arise naturally 



^We explicitly state open the possibility, since the limited resolution of the telescope does not allow us to 
go beyond this possibility. 
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from our model because we are modeling the wind self-consistently with complete treatment 
of radiative processes in all wavelengths. In our models these UV features as well as the 
long wavelength absorption lines in the X-ray band, are produced by the low ionization pa- 
rameter part of the flow (e.g. log.^ = 0.12,-1.13 from Figure [HI lower right panel). This 
demonstrates that the same absorber can produce X-ray and UV lines with similar velocities. 
Figure [12] shows the kinematic relationship between the X-ray and the UV absorbers. The 
model has been shifted (up) for clarity. Here we show the O vill A18.969 and the doublet 
O VI AA1032, 1038 (with respect to the shorter wavelength). For comparison we plot the 
histogram data of the 900 ks of NGC 3783. We see that our model is capable of simulta- 
neously producing X-ray and UV absorption lines with similar velocities around ~ 1000 km 
s~^, as it has been suggested from UV data (iGabel et al.l 120031 ). 



5.3. Present single wind vs Multiphase wind scenario 



However, we want to highlight the most important differences between our model, con- 
sisting of a single wind, governed by the laws of radiative acceleration and the multiphase 
wind; subjected to pressure equilibrium. 

The first difference is that in our model there is a clear correlation between ionization 
state of the ions and velocity {^ — v), and C, and number of particle uh (^ — uh). At the same 
time because of the dependence of v and uh on the spatial distance r, a dependence of ^ on 
r. This is di f ferent from the multicomponent in pressure equilibrium model, suggested by 
Netzer et al.l ( l2003l ). In that model three components (in ionization) may coexist in the same 
volume of space, and lying on the stable regions of the nearly vertical part of the thermal 
stability curve {logT vs log[U/T] see Figure 12 in that paper). Our absorbers are not embed- 
ded in an external medium, and they cannot coexist at ex actly the sarne locat ion but follow 
the physical laws of radiative acceleration. In that context, iGabel et al.l (120051 ) also find that 
the UV absorbers of NGC 3783, may share some properties of the multiphase thermal wind. 
The UV kinematic components lb, 2 and 3 could occupy the low-temperature base of the 
region of the thermal stability curve where a range of temperatures can coexist at pressure 
equilibrium. There is one weak point in the picture of inhomogeneities coexisting at pressure 
equilibrium. The low-ionization high- velocity UV absorber does not fit there, due to a factor 
of 10 larger pressure than the other component. If embedded into a more ionized- hotter 
material, it will eventually evaporate, unless there exist an additional confining mechanism. 
This extra confining could be provided by magnetic pressure, requir ing moderate magnetic 
fields (B ^ 10~^ G), as i s predicted by some dynamical models ( lEmmering et al.l Il992 



de Kool &: Begelmanlll995[ ). In our model, there is no need for confining mechanism, and if 
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any the more ionized material is closer to the source shielding the intense continuum radi- 
ation preventing its evaporation. We conclude by stating that our model is one competitor 
more, it cannot be ruled out compared with pressure equilibrium model. A more detailed 
comparison between models is beyond the scope of the present work, that could take place 
somewhere in the future. 



6. Two outflows 

So far we have been concerned with the line velocity shifts as measured with respect to 
the points of maximum absorption of the lines, but little have been said about the asymmetry 
of the troughs. To quantify the line asymmetries predicted by the theoretical models we have 
tabulated the fraction of the terminal velocity at which the absorption is maximum, to be 
compared with the ratio {v/vi), where vi is the position of the blue edge of each line, used 



to qu antify the asymmetry in the spectrum of NGC 3783, see Table 3 of iRamfrez et al. 



(120051 ). We characterize a theoretical line as "red" if v{Tmax)/'Voo > 0.5 and as "blue" if 
v{Tmax) /voo < 0.5, while for observed lines they are grouped with respect to v/vi > 0.5 and 
v/vi < 0.5 (see Ramirez et al. 2005). In Table |3]we present the observed troughs. In the first 
and second columns are the identifications and the wavelengths of the lines, in the third we 
give the fraction v( Tm.nT) /vrx^ for model A , and in the fourth column the classification from 



observed lines, as in lRamirez et al.l (|2005[ ). i.e., v/vi. We classify every trough as either (R) if 
it is red or (B) if it is blue. There are clear discrepancies between model A and observation. 
In order to improve on the agreement with observations we found necessary to create a 
model composed by two outflows, which differ by their log^o^ the initial exposition to the 
source. In table H] we can see the parameters of the two outflows separately, which composed 
model C. The first flow with log^o = 3, which we call HIF (stand from high ionization 
flow), is able to create the majority of the high ionization lines (S xvi A4.729, S XV A5.039, 
Fe XXIII A8.303), with extended blue wings as the flow moves outwards, and acquires more 
velocity, in agreement with the observations. The second, which we call LIE (stand from low 
ionization flow), is characterized by log^o = 2, is able to create the extended blue wing of 
the oxygen lines. This model has the fundamental ingredients that explain the asymmetry 
observed in the lines of NGC 3783 and at the same time, gives us the bulk velocity of the 
flow. 

In the third column of Table E] we present the characterization of the lines coming from 
the two-outflows model. The HIF component, is the major contributor to the formation of 
the high ionization lines, i.e. S xvi, Fe xxi, Fe xxil, Fe xxill, because their optical depths 
peak close to log^o(HIF), and as the velocity of the flow increases their ionization fraction 
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decrease, forming their blue wings. Under this picture, only a few lines are characterized 
as blue, while the rest are red. Some intermediate-to-low ionization ions, have appreciable 
fractions in the LIF component, which broadens the lines and yields a red apparency to the 
composed profile. This effect is illustrated by the fourth column of Table O 

The composed model reproduces the asymmetry of the lines. The present uncertainties 
in the computation of v{Tmax) are of the order of 60 km s~^ (or Av{Tmax)/voo = 0.055 with 
foo = 1100 km s~^). So the lines classified as pure red with vijmax) > 0.55 and pure 
blue with v{Tmax) hoo ^ 0.45 are more reliable than those with 0.45 < vijmax) hoo ^ 0.55. 

It is interesting to analyze why the O Vll and O vill lines are blue, while the lines 
Mg XI A7.473 and Fe xvil A15.015, remain red. The optical depths of Mg xi A7.473, and 
Fe XVII A15.015 are below 10~^ at velocities greater than 1100 km s~^, so the LIF component 
does not contribute to their broadening. Another interesting example is the difference in 
classification between the lines Mg xi A 7. 473 and Mg xi A9.169. The reason is that while 
for Mg XI A7.473 the oscillator strength is ~ 5 x 10~^, the f-value for the line Mg xi A9.169 
is almost one order of magnitude higher (~ 7 x 10^^), which makes the optical depth of the 
latter significant up to velocities beyond ~ 1600 km s~^. 

Figure [13] compares the modeled velocities of the lines from the two-outflows model with 
observations. Robust regression using the M estimator gives a slope of 0.63±0.22 (solid line 
in the figure), and the residual standard error (rms) is 16.6 for 25 degrees of freedom (dof). 

A comparison between the predicted trough profiles and those observed in NGC 3783 
is shown in Figure [TH Here we plot the theoretical spectrum given by the outflows HIF and 
LIF, in the 18.5 — 19.3 A. Here we have used a powerlaw as continuum with spectral index 
of F = 1.77, suitable for this AGN, set an extra absorber the wabs @ model in XSPEC, 
and the composed fluxes generated by the HIF and the LIF components. To compose the 
O VIII A18.969 line profile we have summed the fluxes from both flows using the equation 
f lT2|) (where m = 22). With this approximation th e profile predic t ed by the two-outflows 



model is in excellent agreement with the observation iRamirez et al.l ( 120051 ). The red wing of 
the troughs, in the range 0-1000 km s~^, is formed by the HIF, while the blue wing, in the 
range ~ 1000-2500 km s~^, is formed by the LIF. These values are in agreement with the 
values measured for this lines in the observed spectrum. 

It is the first time that such theoretical work is performed to explain the asymmetry 
seen in absorption in the X-ray spectrum of AGNs. Other works have been able to explain 
the blue wings of UV lines, like the well-studied C IV A 1549, for example from radiation- 



^Nwabs = 0.1 X 1022 Cni-2. 
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driven disk- wind models ( lProgall2003l ). So we cannot conclude that the model presented here 
is unique for the description of the asymmetry seen in X-ray lines of AGNs, until detailed 
comparisons are made with suc h models, including X-ray line profile produced by a wind 



from a Keplerian acc retion disk (iKnigge et al.lll995l : IShlosman et al.lll996l : iProga et al.ll2000 



Drew fc Progal 120001 ) ■ Further work is necessary to be performed for such comparison and 



we plan to do that in a near future. 



7. Summary 



We have computed photoionized wind-flow models for the X-ray spectrum of NGC 3783. 
We studied singly continously absorbing models as well as multiple optically thin compo- 
nents linked through an analytic wind velocity law. It is found that the singly continuously 
absorbing model yields gas column densities and optical depths too high, unless one adopts 
very low (~ 10^~^ cm~^) densities and metal abundances (~ 0.01 solar). On the other hand, 
the multicomponent model is able to reproduce observations very well. For this model we 
compute ionization properties of the material using a velocity law compatible with a radiative 
wind. Our model is consistent with logrio ~ 11.35 [cm~^], a launching radii of logro ~ 15 
[cm], and a terminal velocities of ~ 1500 km s~^, which yields a mass loss rate of the 
order of Mout ~ 1 ^Q/yr (assumin g a volumic factor = 0.1). If we assume an ionizing 
luminosity of L ~ 2 x 10'^^ ergs s~^ (jPeterson et al.ll2004l ). and accretion efficiency of rj = 0.1, 
th e Eddington mass accr e tion m ass is Medd ~ 0.01 M(T^/yr. This is c onsistent with the result 
of ICrenshaw fc Kraemeil J2007h for NGC 4151, and iRamfre j J2008h for A PM 08279+5255 
of M^ut/Medd > 10. However, it is different from the supposition made by iGoncalves et al. 



( l2006l ). of Mout/Medd < 1, using their photoionization code TITAN, for computing the single 
medium in pressure equilibrium. 

Finally, the asymmetry seen in the lines of the X-ray spectrum of NGC 3783 required a 
model with two outfiows. One flow with a launching ionization parameter of log^o ~ 3 and a 
column density of Nh = 10^^ cm~^, which recreates the red wings of the low ionization lines 
from Ne IX to O vil, and the blue wings of the high ionization lines from Fe xxill to Si xiv. 
A second flow is necessary to to create the blue wing of the oxygen lines, which exhibit a 
blue character, the theoretical fltting required log^o ~ 2, terminal velocities of around ~ 
2200 km s~^, and a column density of A'^^^ 



10^1 cm-2. 



Our calculations also predict a relationship between the UV and X-ray bands, as models 
adjusted to flt the X-ray spectrum naturally predict UV line s like the Lyman serie, and the 
O VI, N V and C iv doublets, in apparent concordance with ICostantinil (120101 ). 
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Table 1. Composition and parameter values of the kinematics model A. 



Element 


Relative Abundance 


H ■• 


1.0 


He ■•■ 


0.1 


C •■• 


0.3540E-03 


N ••• 


0.9330E-04 


••• 


0.7410E-03 


Ne •■• 


0.1200E-03 


Mg •■• 


0.3800E-04 


Si ■•■ 


0.3550E-04 


S ••■ 


0.2140E-04 


Ar •■• 


0.3310E-05 


Ca ••• 


0.2290E-05 


Fe •■• 


0.3160E-04 



'^In this model wq = 0.4, Voo = 
900 km s^^, log;^Qro = 15.75 [cm], 
Vturb = 200 km s~^ and Nh = 5 x 
10^0 cm-\ 
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Table 2. Composition and parameter values of the kinematics model B. 



Element 


Relative Abundance 


H ■• 


1.0 


He ■•■ 


0.1 


C •■• 


0.3540E-03 


N ••• 


0.9330E-04 


••• 


0.7410E-03 


Ne •■• 


0.1200E-03 


Mg •■• 


0.3800E-04 


Si ■•■ 


0.3550E-04 


S ••■ 


0.2140E-04 


Ar •■• 


0.3310E-05 


Ca ••• 


0.2290E-05 


Fe •■• 


0.3160E-04 



^In this model wq = 0.4, Voo = 
1100 km s^"*^, logiQ Vq = 15.75 [cm], 
Vturb = 200 km s~^ and Nh = 5 x 
10^0 cm-\ 
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Table 3. Asymmetry comparison between model A and observation 



Ion Line (A) v{Tmax) /QOO km s'~'^ Ramirez et al. f2005 ) " 



S XVI 


4.729 





55 


(R) 


0.33 ± 0.1 (B) 


b XV 


5.039 





67 


(R) 


0.23 ± 0.1 (B) 


Si XIII 


5.681 





70 


(R) 


0.35 ± 0.1 (B) 


Si XIV 


6.182 





67 


(R) 


0.23 ± 0.08 (B) 


Si XIII 


6.648 





70 


(R) 


0.43 ± 0.074 (B) 


Mg XII 


7.106 





70 


(R) 


0.32 ± 0.070 (B) 


Mg XI 


7.473 





76 


(R) 


0.74 ± 0.067 (R) 


Fe XXIII 


8.303 





55 


(R) 


0.037 ± 0.060 (B) 


Mg XII 


8.421 





70 


(R) 


0.32 ± 0.059 (B) 


Mg XI 


9.169 





76 


(R) 


0.31 ± 0.054 (B) 


AT r 

Me X 


9.708 





76 


(R) 


0.45 ± 0.052 (B) 


Ne X 


10.240 





76 


(R) 


0.17 ± 0.048 (B) 


AT , 

Ne IX 


11.547 





82 


(R) 


0.69 ± 0.043 (R) 


Fe XXII 


11.770 





55 


IT) \ 

(R) 


O.zz ± 0.042 (Bj 


Ne X 


12.134 





76 


(R) 


0.29 ± 0.041 (B) 


Fe XXI 


12.284 





62 


(R) 


0.42 ± 0.040 (B) 


Fe XVIII 


14.373 





70 


(R) 


0.27 ± 0.035 (B) 


Fe XVIII 


14.534 





70 


(R) 


0.51 ± 0.034 (R) 


VIII 


14.832 





82 


(R) 


0.21 ± 0.033 (B) 


Fe XVII 


15.015 





70 


(R) 


0.76 ± 0.033 (R) 


VIII 


15.188 





82 


(R) 


0.38 ± 0.033 (B) 


VIII 


16.006 





82 


(R) 


0.26 ± 0.030 (B) 


VII 


17.200 





90 


(R) 


0.29 ± 0.029 (B) 


VII 


17.396 





90 


(R) 


0.30 ± 0.029 (B) 


VII 


17.768 





90 


(R) 


0.32 ± 0.028 (B) 


VII 


18.627 





90 


(R) 


0.32 ± 0.027(B) 


VIII 


18.969 





82 


(R) 


0.26 ± 0.026 (B) 



f/fi, see definition in the text. 
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Table 4. Parameters of the two outflows model (Model C) 



Outflow HIF Outflow LIF 



logCo^S.O logeo = 2.0 

= 1100 km s-^ = 2200 km s"^ 

Nh = 10^2 cm-2 Nh = 10^1 cm-2 

Wo = 0.4 Wo = 0.2 

Solar composition Solar composition 



Table 5. Asymmetry comparison between model C and observation 



Ion Line (A) t;(t^„^[HIF])/1100 km s"^ t;(t^„^[HIF+LIF])/2200 km s"^ 



S XVI 


4.729 


0.47 


(B) 


0.47 


(B) 


S XV 


5.039 


0.55 


(R) 


0.55 


(R) 


Si XIII 


5.681 


0.62 


(R) 


0.62 


(R) 


Si XIV 


6.182 


0.55 


(R) 


0.55 


(R) 


Si XIII 


6.648 


0.62 


(R) 


0.62 


(R) 


Mg XII 


7.106 


0.55 


(R) 


0.55 


(R) 


Mg XI 


7.473 


0.74 


(R) 


0.74 


(R) 


Fe XXIII 


8.303 


0.47 


(B) 


0.47 


(B) 


Mg XII 


8.421 


0.55 


(R) 


0.55 


(R) 


Mg XI 


9.169 


0.74 


(R) 


0.37 




Ne X 


9.708 


0.74 


(R) 


0.74 


(R) 


Ne X 


10.240 


0.74 


(R) 


0.37 




Ne IX 


11.547 


0.74 


(R) 


0.37 




Fe XXII 


11.770 


0.47 


(B) 


0.47 


(B) 


Ne X 


12.134 


0.74 


(R) 


0.37 


(BY 


Fe XXI 


12.284 


0.47 


(B) 


0.47 


(B) 


Fe XVIII 
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With the addition of the LIF, we mark those line whose symmetry change. 
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Fig. 1. — Statistical significance of tlie sliifting of the X-ray lines found in the spectrum of 
NGC 3783. Each box is made of five-number summaries: the smallest observation (sample 
minimum, which are the extreme thin bars at the left of each box), a lower quartile (Ql, 
which is the left thick border of each box), median (Q2, black thick vertical line), upper 
quartile (Q3, the right thick border of each box) and largest observation (sample maximun, 
which are the extreme thin bars at the right of each box). 
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Fig. 2. — Optical depth (r) static case vs outflow model for the line Ne x A12.134 (see text). 
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Fig. 3. — Profile of tlie line Ne x A12.134 formed by taking F\{v) = Fc x exp[— TA(t')] as 
function of the velocity. The rectangle profile is made by a cloud with solar abundances 
(open circles). The profile with an extended red wing and a sharp blue wing is made by a 
cloud with 1% solar abundances (open square). 
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(x)iw nOTx 

Fig. 4. — Variation of the main kinematic and ionization variables for model A (see text). 
Upper left: Velocity law w{x) for the model. Upper right: The ionizing radiative flux, 
decaying with oc r~^ (dashed line) and with oc r~^~^ (dotted line). Lower left: Variation of 
the density as function of x. Lower right: Variation of the ionization parameter (^) with the 
normalized velocity. 
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Fig. 5. — Measured vs model p redicted velocities for the group of unblended li i ies gi ven in 
Table 2 of iRamirez et al.l ( 20051 ). The measurements are taken from lKaspi et al.l ( 20021 ) . The 
solid line is the best line after a linear regression for model A. 
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Fig. 6. — Variation of the main kinematic and ionization variables for model B (see text). 
Upper left: Velocity law w{x) for the model. Upper right: The ionizing radiative flux, 
decaying with oc r~^ (dashed line) and with oc r~^~^ (dotted line). Lower left: Variation of 
the density as function of x. Lower right: Variation of the ionization parameter (^) with the 
normalized velocity. 




Fig. 7. — Measured vs model p redicted velocities for the group of unblended li i ies gi ven in 
Table 2 of iRamirez et al.l ( 20051 ). The measurements are taken from lKaspi et al.l ( 20021 ) . The 
solid line is the best line after a linear regression for model B. 
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Fig. 8. — Relationship between the ionization parameter and velocity. The solid line is 
the theoretical predi ction for model A. The open squares with error bars are points with 
velocities taken from lKaspi et al.l (|2002[ ). for the group of lines compared in Figure O 
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Fig. 9. — Relationship between the ionization parameter and velocity. The solid line is 
the theoretical predi ction for model B. The open squares with error bars are points with 
velocities taken from lKaspi et al.l (|2002[ ). for the group of lines compared in Figure [71 
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Fig. 10. — Predicted UV spectrum by model B (in the range 700-1700 A). These absorption 
lines are formed by the faster (or at lower ionization parameters) clouds of model B, i.e., 
log^ — 0.12, —0.65. Most of these hues are found in real UV spectrum of AGN. 
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Fig. 12. — Possible kinematical relationship between the UV and the X-ray absorbers. In 
the figure are plotted the theoretical X-ray and UV hnes O Vlll A18.969 and O VI A1032, 
dashed and short-dotted lines respectively. Also plotted the profile of the line O Vlll A18.969 
taken from the 900 ks spectrum of NGC 3783. All these hnes are sharing velocities ~ 1000 
km s~^. 
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Fig. 13. — Measured vs mo c iel pr edicted velocities for the group of ur iblended lines give n 
in Table 2 of iRamirez et al.l ( l2005l ). The measurements are taken from iKaspi et al.l ( 12002| ). 
The solid line is the best line after a linear regression for model C. 
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Fig. 14. — A graphical representation of the flows HIF and LIE with the line profiles seen 
in NGC 3783 (histogram). Solid line is the composed model (model C) by the two flows. 
It can be seen that it was necessary to include the second low ionization flow in order to 
reproduce the asymmetry of these lines. Otherwise the high ionization flow would form a 
more extended red wing. 



